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Abstract 



An enhanced grid system for the Orbiter was built by integrating CAD definitions from 
several sources and then generating the surface and volume grids. The new grid system 
contains geometric components not modelled previously plus significant enhancements on 
geometry that has been modelled in the old grid system. The new Orbiter grids were then 
integrated with new grids for the rest of the launch vehicle. 

Enhancements were made to the hyperbolic grid generator HYP GEN and new tools 
for grid projection, manipulation and modification, Cartesian box grid and far field grid 
generation and post-processing of flow solver data were developed. 



Introduction 

This project forms part of the long term computational effort to simulate the time 
dependent flow over the integrated Space Shuttle vehicle during its ascent mode for various 
nominal and abort flight conditions. One of the important parts of this project is to obtain 
better comparisons between numerical simulations and flight data. In the simulations 
performed so far, the geometry has been simplified in various ways to reduce the complexity 
so that useful results can be obtained in a reasonable time frame due to computer resources 
limitations. With the ICEM CAD/surface grid generation software, the surface definitions 
of different components of the Space Shuttle are now more accurately represented. In this 
project, surface and volume grids with more refined details are generated for the Orbiter. 
The complex process of inter-grid connectivity is then performed using the PEGASUS code. 
The robustness of the new grids in the Orbiter base region is tested by first combining with 
the old Orbiter fuselage and wing grid and performing a flow simulation for an Orbiter 
alone case. The enhanced new Orbiter grid system is then integrated with the enhanced 
External Tank, Solid Rocket Boosters and attach hardware grid systems developed at 
NASA-Johnson to form the enhanced Space Shuttle Launch Vehicle (SSLV) grid system. 

In the past year, the hyperbolic grid generation code HYPGEN has been modified 
significantly to achieve a more modular structure. The algorithm has also been improved 
to provide robustness for a wide range of different geometries [1]. This project includes 
further research in developing even more robust algorithms. The ultimate goal is for the 
code to be able to intelligently select as many input parameters related to grid quality as 
possible. Thus, the user input will be reduced to a minimum. 

Various useful tools are developed to improve and speed up the process of grid genera- 
tion and manipulation and flow solution data post-processing. These include improvements 
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made to the grid projection tool PROGRD used to project subsets of an active grid on to 
subsets of reference grids. This is an essential step in the Chimera overset-grid approach 
for inter-grid connectivity at the body surface. A grid editor GRIDED is developed that 
includes many common operations performed on a single block structured grid into a uni- 
fied program. Another program called BOXGR is developed to generate the Cartesian box 
grids and far field grid used in the Chimera overset grid approach (see section on BOXGR 
below). A post-processing tool OVERPOST is developed to convert the output data files 
generated by the flow solver OVERFLOW [2] into a format suitable for plotting using 
XYPLOT. 

Algorithm and Code Development 

Hyperbolic Grid Generator - HYPGEN 

Further optimization was performed on the code including the development and im- 
plementation of more efficient block tridiagonal solvers resulting in a 10% speed up of 
the code. Slight algorithm improvements were made for the sharp convex corner logic 
and floating boundary conditions at singular axis points. Version 1.2 of the code was 
released in March, 1992 to users in the RF Division at Ames and to the Shuttle group 
at NASA Johnson. In 1992, the code was requested by and distributed to the follow- 
ing groups: Rockwell International, Nielson Engineering, Computational Fluid Dynamics 
Applications, Penn State University, The Boeing Company, Wright Lab., Boeing Defense 
and Space Group, Office of Aircraft Compatibility (Dept. AF), Douglas Aircraft, NASA 
Langley, Bell Helicopter, B-2 Division (Northrop Corp.). A detailed description of the 
algorithm for HYPGEN is given Appendix A. 

Grid Projection Tool - PROGRD 

This program reads in an active grid and performs projections of subsets on the body 
surface (L=l shell) of the active grid on to subsets on L=1 of one or more reference grids. 
Projections can be performed along surface normals or along any of x, y or z directions. 

PROGRD is required where solid surfaces from two grids overlap. This occurs not 
only in the overlapped regions between a collar grid and its parent grids but also in complex 
overlapped grid regions such as the backend of the orbiter. There are two main areas where 
PROGRD is required: 

(1) The outer boundary of one grid has to lie on the bilinear surface of the neighboring 
grids for PEGASUS to find the right interpolations, e.g., the outer boundary of the 
collar grid has to be projected on to the L=1 surfaces of the parent grids. 

(2) Hole boundary points of one grid has to lie on the bilinear surface of a neighboring 
grid for PEGASUS to find the right interpolations, e.g., the hole boundary points of 
a parent grid have to lie on the L=1 surface of the collar grid. Version 3.2 or later of 
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PROGRD has the added ability to look for fringe orphan points with an iblank value 
of 101 (compatible with latest PEGASUS 4.1 +). 

The later versions of PROGRD (3.0 or later) allows the recording all projections that 
have to be performed for the active grid under a single input file. Not all subsets recorded in 
the input file axe projected. Projections are performed only on the subsets whose reference 
grid file can be found under the specified directory. This allows continuous modifications 
of the active grid without keeping the old reference grids around. Version 3.1 or later also 
accepts a different mode and projection direction for each subset. 

Grid Editing Tool - GRIDED 

This is an effort to include many common operations that axe performed on a grid 
into a single code. Separate programs are in existence to perform each of these tasks but 
it is cumbersome to have to carry all the different programs around. 

The GRIDED program can be used to perform the following common operations on a 
2D/3D single grid file in PLOT3D format. The operations can be performed in sequence 
one after another. After each operation, a warning message is printed if the original grid 
has changed from a right-handed to a left-handed system or vice versa. The grid indices 
are labelled as J, K and L. 

Operations that do not change the grid dimensions include 

(1) Convert to/from formatted/unformatted, planes/whole, 2D/3D data types 

(2) Interchange J and K grid families 

(3) Interchange K and L grid families 

(4) Interchange J and L grid families 

(5) Reverse index direction in J and/or K and/or L 

(6) Scale and translate 

(7) Mirror about x=0, y=0 or z=0 plane 

Operations that change the grid dimensions include 

(1) Extract a subset 

(2) Add any number of extra planes in J, K or L by extrapolation 

(3) Add extra planes in J or K by reflected symmetry 

(4) Add extra planes to 2D grid to form 3D grid for ‘2D’ or ‘axisymmetric’ option in 
OVERFLOW 

(5) Duplicate and mirror reflected symmetry grid to form periodic grid 

(6) Read in a new grid and concatenate to current grid 

Box and Far Field Grid Generation Tool - BOXGR 

In the Chimera overset grid approach for complex configurations, it is typical to grow 
the volume grids only a small distance out from each componenet and then embed the 
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component volume grids in a Cartesian box grid. The Cartesian box grid typically consists 
of an uniformly spaced interior region which completely encloses the component volume 
grids, and outer layers where the grid spacing is stretched from the interior value to some 
desired value at the outer boundary. A tool called BOXGR was developed to generate 
box grids with arbitrary dimensions and grid spacings in the interior and outer regions. 
The same tool also has the capability to generate the inner surface of the fax field grid 
where at least one point overlap is maintained with the outer boundary of the box grid. 
An ellipsoidal topology, built from smoothing the corners of the box grid, is used for the 
inner surface of the fax field grid. 

Post-processing Tool for OVERFLOW - OVERPOST 

A post-processing tool named OVERPOST was developed for the Chimera overset- 
grid Navier-Stokes solver OVERFLOW. The function of OVERPOST is to extract appro- 
priate data from the output files from OVERFLOW for plotting purposes using XYPLOT. 
Processed files axe generated for histories of conservative variables residuals, forces and mo- 
ments, minimums of pressure and density, turbulence model equations residuals and species 
residuals. A preference file can be initialized such that the user does not have to enter the 
same parameters each time OVERPOST is used. 

Applications 

The construction of the new enhanced Orbiter grid system is part of a team effort 
with NASA Johnson to enhance the grid system for the complete Space Shuttle Launch 
Vehicle. The first part of the work for this project is to generate surface grids for the 
Orbiter base region. Geometric components that were absent from the old grid system are 
now modelled. These include the body flap, the Space Shuttle Main Engines (SSMEs), 
the OMS extension and the RCS jet block. The surface geometry was obtained from 
Rockwell and ICEMCFD was used to manipulate the geometry and generate the surface 
grids. One region where geometry simplification was needed was in the RCS jet block 
where the CAD surface was modified to join smoothly with the OMS extension. Volume 
grids were generated using HYPGEN, and the inter-grid connectivity was performed using 
PEGASUS by combining the new base grids with the old Orbiter fuselage and wing grid 
resulting in an intermediate grid system consisting of 18 grids. The robustness of the new 
base grids was tested by running a flow solution using OVERFLOW with the intermediate 
18-grid system (see Figure 1). Since the solution was fairly smooth in the base region, it 
was concluded that the topologies used for the new base grids were adequate to capture 
the flow in the region. 

The next step was to generate enhanced grids for the Orbiter fuselage, wing and 
elevons. In the old grid system, the Orbiter wing and elevons were generated with some 
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coarse approximations in places and no CAD definition was in existence at that time. 
Part of the effort in this project was spent on the construction of CAD surfaces for part 
of the wing from the Masters Dimensions Specifications [3]. During the course of the 
work, CAD surfaces for the wing tip, the elevons and parts of the wing were obtained 
from Kennedy Space Center and Rockwell. These latest CAD definitions were integrated 
together with the existing CAD model to form enhanced geometry definitions for the wing 
and elevons. A blocking strategy was developed for the model which replaced the single 
old Orbiter fuselage/ wing grid with about 11 overlapped grids for the right half: Orbiter 
nose (ONOSE), Orbiter forward fuselage (0RB1), Orbiter back fuselage (0RB2), upper 
wing (WINGU), lower wing (WINGL), wing leading edge (WINGLE), elevons at 10/5 
deflections (ELVS105), middle elevon gap (EGAPM105), in-board elevon gap (EGAPI10), 
Orbiter lower-middle fuselage (ORBLFM), Orbiter lower-back fuselage (ORBLFB) - see 
Figures 2 and 3. Introduction of the ONOSE grid removed the axis singularity at the nose. 
The WINGU and WINGU grids were designed to have grid lines alined with pressure taps 
where flight data is available for comparison with CFD results. The separate elevon and 
elevon gap grids allow different deflection angles to be modelled without changing the grids 
on rest of the Orbiter. A fine ORBLFM grid was used to enclose the complex grids for the 
ET/ Orbiter attach hardware. 

All volume grids were generated using HYPGEN except for the elevon gaps which 
were generated using ICEMCFD with a patched grid topology. The following box grids 
were generated using the BOXGR tool described above: a large box around the entire 
Orbiter with a 12 inch interior grid spacing and stretched to a 48 inch grid spacing at 
the outer boundaries, a small box around the Orbiter base grids with a uniform 5 inch 
grid spacing and a small box around the tip of the vertical tail. Inter-grid connectivity 
was performed using PEGASUS 4.1 for the right half of the Orbiter 30 grids). The 
right half grids were then reflected to form left half grids. Some of the right and left grids 
were combined to form a single grid. Inter-grid connectivity was then performed for the 
complete Orbiter (« 50 grids). 

The final stage of the grid generation was to integrate the new enhanced Orbiter 
grid system with enhanced grids for the rest of the launch vehicle. Another box grid was 
generated to enclose the ET/SRB and the far field grid was also generated using BOXGR 
and HYPGEN. Slices of the SSLV volume grids, box grids and the far field grid on the 
symmetry plane are displayed in Figure 4. A 5-day meeting was held at Ames with the 
team from NASA-Johnson and Lockheed, Houston to integrate the enhanced SSLV grid 
system using PEGASUS 4.1. The final grid system contained 111 grids with 16 million 
points. More details of the development of this grid system are discussed in Appendix B. 
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Concluding Remarks 


Useful tools were developed for grid projection (PROGRD), grid manipulation and 
modification (GRIDED), box and far field grid generation (BOXGR) and post-processing 
of OVERFLOW data files (OVERPOST). These tools have been distributed to users in the 
RF Division at Ames and to the Shuttle group at NASA-Johnson. Further optimization 
was performed on HYPGEN together with some slight algorithm improvements. HYPGEN 
has been distributed to many different groups in industry. 

Enhanced CAD geometry definitions for the Orbiter obtained from several sources 
were combined into a single model. New grids were generated to model components absent 
from the previous grid system and to enhance resolution of existing components. The new 
Orbiter grids were integrated with new grids for the rest of the launch vehicle. The 
complete enhanced integrated SSLV grid system with 111 grids and 16 million grid points 
is one of the largest overset grid system built for complex aerodynamic configurations. 


References 

1. Chan, W.M. and Steger, J.L., “A Generalized Scheme for Three-Dimensional Hy- 
perbolic Grid Generation”, AIAA Paper 91-1588. Proceedings of the AIAA 10th CFD 
conference, Honolulu, Hawaii, 1991. 

2. Buning, P.G., Chan, W.M., Renze, K.J., Sondak, D., Chiu, I.T. and Slotnick, J.P., 
OVERFLOW User’s Manual, Version 1.6, NASA Ames Research Center, Moffett Field, 
CA, 1991. 

3. “Space Shuttle Masters Dimension Specifications”, Rockwell International MD- 
V70, Rockwell Space Systems Division, Downey, CA, Revised 1991. 


6 




ORIGINAL PAGE 
COLOR PHOTOGRAPH 


Figure 1 




SURFACE GRIDS FOR NEW ORBITER GRID SYSTEM 



ORIGINAL PAGE 

COLOR PHOTOGRAPH 


Figure 2 


SURFACE GRIDS FOR NEW ORBITER GRID SYSTEM 



ORIGINAL PAGE 
COLOR PHOTOGRAPH 


Figure 3 






APPENDIX A 



Enhancements of a Three-Dimensional 
Hyperbolic Grid Generation Scheme* 

William M. Chan 

MCAT Institute 

NAS A- Ames Research Center 

Mail Stop T045-2 

Moffett Field , California 94035 

and 

Joseph L. Steger** 

Department of Mechanical Aeronautical and Materials Engineering 
University of California at Davis 
Davis , California 95616 


Transmitted by Joe F. Thompson 


ABSTRACT 

A hyperbolic grid generation scheme formulated from grid orthogonality and cell 
volume specification has been significantly enhanced so that high quality three-dimen- 
sional grids can be obtained for a wide variety of geometries. While the speed of the 
scheme remains one to two orders of magnitude faster than typical elliptic grid 
generation methods, the robustness of the scheme has been greatlv improved over 
previous applications of the three-dimensional hyperbolic grid generation procedure. 
Enhancements included the use of spatiallv-variable smoothing coefficient, metric 
correction procedures, local treatment of severe convex comers, and new extrapola- 
tion treatments of floating and axis boundaries. The versatility of the new hyperbolic 
grid generation scheme is demonstrated by three-dimensional grids generated for 
external components of the integrated Space Shuttle vehicle and the SOFIA tele- 
scope. 


•The authors are grateful to Dr. Pieter Buning and Dr. Yehia Rizk for some very useful 
discussions. A special thanks is owed to Mr. Chris Atwood for providing the surface grid for the 
SOFIA telescope. 

••Joseph Steger died on May 1, 1992; Computation Fluid Dynamics has lost a giant. 
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1. INTRODUCTION ‘ • 

One of the most popular approaches for generating structured grids is bv 
the solution of a set of partial differential equations. The governing equations 
can be classified into three types: elliptic, parabolic and hyperbolic. The most 
widely used grid generation methods require the solution of a set of elliptic 
equations [4, 14-16]; however, parabolic and hyperbolic equations have also 
been successfully employed [3, 6, 8, 10, 12, 13] and are advantageous for 
certain applications. 

Since the solution of the elliptic equations satisfies the maximum principle, 
the grids generated are typically smooth. Moreover, the formulation of the 
elliptic equations allows exact specifications of all boundary point locations. 
However, grid orthogonality cannot be maintained with conventional elliptic 
grid generation methods and boundary surface refinement can be difficult to 
impose. The user input required to set up boundary distributions can also be 
time consuming. Hyperbolic grid generation schemes produce nearly orthog- 
onal grids, have excellent clustering control, and can be generated in one to 
two orders of magnitude less computer time than elliptic methods. However, 
hyperbolic grid generation methods are less robust, tend to propagate input 
discontinuities, and the outer boundary location cannot be precisely specified. 
Hence, they are usually restricted to the generation of grids for external flows 
or for chimera overset-grid schemes [2], where the exact location of the outer 
boundary is not constrained. This latter application has become more impor- 
tant, however, and this, coupled with the efficiency and often superior grid 
quality obtained with hyperbolic grid generation schemes, motivates investi- 
gation into improving their robustness. 

In the last few years, the hyperbolic grid generation algorithm described in 
[13] has evolved to include some significant enhancements. A wider range of 
boundary conditions can now be treated and some feedback features have 
been added so that dissipative terms (which give the equations a somewhat 
parabolic nature and smoothness) are now automatically adjusted depending 
on grid evolution or character. The resulting hyperbolic grid generation 
scheme is significantly more robust, produces higher quality grids, and can 
treat a wider variety of topologies. 

The governing equations for three-dimensional hyperbolic grid generation 
are presented in Section 2. The numerical marching scheme employed to 
solve these equations is described in Section 3. Three factors important in 
controlling gnd quality, the boundary conditions, cell volume specification 
and added smoothing are discussed in Sections 4, 5 and 6, respectively. 
Discontinuities due to body shape or initial grid point distribution can present 
special problems to a hyperbolic solver that dissipation alone cannot satisfac- 
torily cure. Comer points are especially difficult. A metric correction proce- 
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dure which is essential in providing smooth grids at comers with uneven grid 
spacings is described in Section 7. Extra robustness at sharp convex comers 
can be achieved by switching from solving the hyperbolic equations to some 
other equations at the convex comer points. These schemes are described in 
Section 8. Two grid quality checks built-in to the grid generation code are 
described briefly in Section 9. In Section 10, examples from grids produced 
with the new hyperbolic grid generator for various external components of 
the Space Shuttle launch vehicle and the SOFIA telescope are presented. 
Finally, conclusions are given in Section 11. 


2. GOVERNING EQUATIONS 

Generalized coordinates f(x, y, z) y 7?(x, y, z\ £(x, y, z) are sought 
where the body surface is chosen to coincide with £(x, y, z) - 0 and the 
surface distributions of £ — const and 7} — const are user-specified. With 
external aerodynamic applications in mind, the location of the outer boundary 
£(x, y„ z) = is not specified. The governing equations are derived from 
orthogonality relations between f and £, between Tj and f , and a cell volume 
or finite Jacobian J constraint [13]: 

x ( x i + yi i J( + s f=f"°* ( 21a ) 

xr, x c + yvyi + z n z i = 0 ' ( 2 - lb ) 


x ( yn z i + *c y f z v + x v l Ji z ( ~ x ( yi z v ~ x * l JC~i ~ x c l u z t = AV - ( 2 - lc ) 

or, with r defined a s (x, y, z) T 


r ( -r. = 0, 


(2.2a) 


= 0, 


(2.2b) 


<?(*■ y. =) 

<*(f.i7.0 


= = av. 


(2.2c) 


Equation (2.1) comprise a system of nonlinear partial differential equa- 
tions in which x, y, and z are specified as initial data at £ = 0. Local 
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linearization of (2.1) about the* slate 0 results in the system of grid generation 
equations 


A o( F “ F o)f + B 0 {r- r 0 )^ + C„(r - r 0 ) f =/ (2.3) 


with 


X C J/f : f ” 

0 0 0 
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(2.4d) 


We can rewrite (2.3) as 


A 0 F ( + B n r n + C t) ? ( + e (2.5) 

with e = (0,0, AV + 2AV 0 ) r . Now C 0 “‘ exists unless (AV„) 0, (2.5) can 

be rewritten as 


c <7' A »'f + C.7 'BnS + 7 i “ Co 


( 2 . 6 ) 
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Although the algebraic verification is * not trivial, C^A,, and Cq ] B 0 are 
symmetric matrices [13]. This implies that the linearized system (2.6) is 
hyperbolic and can be marched with f serving as the "time-like” direction. 

3. NUMERICAL MARCHING SCHEME 

The system of grid generation equations given by (2.5) are solved with a 
non-iterative implicit finite difference scheme. An unconditionally stable 
implicit scheme has the advantage that the marching step size in f ran be 
arbitrarily selected based only on considerations of accurately generating the 
grid. Linearization is performed about the previous marching step in 

Let = 1 such that f « j — 1, t) = k — L and £ = /-!. 

Central spatial differencing of (2.5) in f and rj with two-point backward 
implicit differencing in £ leads to 

Al&f ( r /+ l * r /) + r l+ \ ~~ r \ ) ~ Cfi; r l- ] “ g/- [ * (3-1) 

where 

0 
0 

and 





Throughout, only those indices that change are indicated; thus, r /4> , =» 
and r^ x r^ x k h etc. 

Multiplying through by Cf 1 and approximately factoring gives 

(/ + cr , B l a,)(/ + cr*A l « < )(r>, l (3.2) 

where I is the identity matrix. The problem is now reduced to solving a 
sequence of block tridiagonal systems. 



186 


WILLIAM M. CHAN AND JOSEPH L. STECER 

Since all f and ^derivatives are approximated bv central differencing 
numerical dissipation terms are added in these directions. For simplicity only 
second differences are used which are explicitly and implicidy included in the 
basic algorithm as 


[/ + Cr'B,8 v - *„(AV),][/ + Cr l A,8 ( - * lf (AV) f ] 

“ F ') = C f‘£*. - K*( AV )f + *„(AV),]r f , (3.3) 


where, for example. 


(AV) T? r — r Jb+1 - 2 r k + r k _ l% 

and with e if = 2e, f and <r iTJ - 2f^. Additional smoothing and implicitness 
pj are achieved by differencing V^F= F as 1 - FJ = (1 -f- - q f 

with 8 0. This differencing in £ is incorporated into (3.3) as 

I' + Ci + Wb,*, -*„(&),] 

xll+(l + e f )Cr'A,S ( -e if (AV) ( ] 

*( F ‘+ 1 “ *0 = C r‘g/ + i - [*, f ( AV )f + £,„(AV),]r;. (3.4) 

The values of 8 { and 0„ are kept at zero unless the body contains concave 
profiles in £ or 77, in which case values of the appropriate 6 of 1-4 are 
effective in preventing grid lines from crossing. 

The coefficient matrices A h B, and C, contain derivatives in 77 and £. 
The derivatives in £ and 77 are obtained by central differencing while the 
derivatives in £ are obtained from (2.1) as a linear combination of f- and 
77-derivatives as follows 


1: 


AV 

Det(C) 


y( z v- l Jv z (' 

x v z l ~ X ( Z V 
x ( Vv - x v!/(j 


= C~‘~ 


(3.5) 
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Det(C) = (y f =„ “ J /,= f ) 2 + (x„; f -x f =„) 2 + (x f y„ - x^) 1 . 

A discussion on a more appropriate way to compute these ^-derivatives is 
given in Section 7. 

4. BOUNDARY CONDITIONS 

Five types of implicit boundary conditions have been implemented in the 
grid generation code at the f and 77 boundaries (except for the axis condition, 
which is only implemented in f). In the following, coordinate increments 
Ar= r j+l — rj are represented by (Ax, A y, A z) T . 

1. Periodicity: All derivatives at the end points in the periodic direction 
are evaluated by 'wrapping around.' A periodic block-tri diagonal solver is 
used for the inversion of the appropriate factor in the left-hand side of (3.4). 

2. Constant Cartesian plane: If a f or 77 boundary is restricted to an 
r = const, y = const or : = const plane, then that value is enforced and the 
other variables are ‘floated/ For example, for an x = const plane at the j = 1 
boundary, x is held constant and y and r are floated using 


Ax | 


' 0 ' 


= 

*!/ 

Us j 

j- 1 

* As j 


3. Symmetry plane: Conventional reflection planes are used to impose 
symmetry about any x = 0 , 1 / - 0 or : = 0 plane and values are updated 
implicitly. For example, to update a reflected plane at j = 1 for a symmetry 
condition about x = 0 corresponding to j = 2 , x reflects odd and y and z 
reflect even as: 


'Ax 1 



f-Ax] 

Ay 


= 

Ay 

k As j 

J-l 


As J 


Since the j — 2 plane in this example may deviate very slightly from x = 0 
away from the body surface due to round-off error, post-processing is done to 
set the x coordinate exactly to zero at j = 2 at each incremental level in /. 



188 


WILLIAM M. CHAN AND JOSEPH L. STEGER 

. 

4. Floating edge : Much as 1 '™ the case of a constant Cartesian plane an 
entire f orjj boundary can be floated using the simple hyperbolic equation 
r (C = 0 or = 0 to update a boundary plane. This is essentially a zeroth- 
order extrapolation of A r from the adjacent interior value, and it often works 
well. Occasionally, however, the floating boundary plane itself may tend to 
"roll in' or ‘lank’ while its neighboring planes in the interior remain smooth. 
Since the interior points are fine, this problem has been remedied by the 
addition of a fictitious line of points on the body surface next to the floating 
edge by linear extrapolation from the interior. The dimension ofthe surface 
grid in a particular direction is temporarily increased by one or two depend- 
ing on whether one or two floating edges are present in that direction. The 
3D grid is then generated over the extended surface grid, after which the 
fictitious plane or planes of points is removed. The addition and removal of 
these fictitious points are carried out internal to the code and is not a concern 
for the user. 

When using the chimera overset-grid scheme, it is particularly desirable to 
have the floating edges splay outwards, i.e., in the direction away from the 
interior of the grid, thus providing better overlap between neighboring grids. 
Although a free floating edge using zeroth-order extrapolation may bend 
inwards or outwards, it is found that using a mixed zeroth- and first-order 
extrapolation scheme tends to bend the edge outwards. For example, at the 
7 = 1 boundary, we have 


'Ax ’ 
Ay 

— 

'Ax 1 
Ay 

+ £ 'X 


[Ax] 

Ay 


(Axl 

Ay 


i A z j 


UsJ 

j~2 

. 

l Ax J 

>2 

II 



where 0 < e ex < 1. Zeroth-order and first-order extrapolation schemes are 
recovered at the two limits of e ej respectively. TvpicaJlv, a value of s ex = 
0.05-0:2 is used. An example using the floating edge condition with e ex =* 0.2 
at the boundaries of a flat plate is shown in Figure 1. 

5. Axis : When the axis logic is used in the j- direction, it is assumed that 
the boundary condition in the ^-direction is either periodic, or that of 
symmetry, or constant planes at both ends. For example, one may generate a 
grid for a complete ellipsoid (periodic in k ) or just one half of it (symmetry or 
constant planes at end points in Jc). 

The treatment of the axis requires special attention in order to produce 
smooth results. Two methods have been developed. The first method re- 
quires the user to adjust certain input parameters and is able to produce 
smooth results for all cases encountered so far. The second method solves the' 
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FiC. 1. Outwurd-spbving edges of a flat plate. 


governing equations on a local coordinate system spanning points to either 
side of the axis. It has the advantage of having fewer parameters to adjust, but 
it is not appropriate for a conical-like axis and will not be described here. 

The first method involves using a mixed zeroth- and first-order extrapola- 
tion and volume scaling. The axis point is updated implicitly bv imposing that 
Ar at the axis is extrapolated by a mixed zeroth- and first-order scheme 
similar to that given by (4.3). The resulting k max predicted values, where 
Jt max is the number of points around the axis, are averaged to produce a 
unique value at the axis. The methods of volume specification described in 
Section 5 below usuallv produce volumes that are too large near the axis. 
Hence, the volumes in the ring of points around the axis are scaled by a 
reduction factor in the range 0.1- 1.0. Typical values of the extrapolation 
factor and volume scaling factor are 0.4 and 0.3, respectively. Figure 2 shows 
the svmmetry plane of an external tank grid with an axis at the pointed nose 
and an axis at the fiat back. Smoothness is maintained at both the front and 
the back regions bv application of the above scheme. 

5. CELL VOLUME SPECIFICATION 

With the hyperbolic grid generation method, one of the means of control- 
ling the grid is by specification of the cell volumes, AV^ * /. Through the cell 
volumes, the extent and clustering of the grid can be modified. Since the cell 
volume at each point must be given, it is clear that the user must devise a 
simple global method for specifying volumes. 

One method is that the specified volume at each point is set equal to the 
computed surface area element times a user-specified arc length. Specifically. 

(5.1) 


ORiGiF*AL PACE !S 
OF POOR QUALITY 
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the Wk' SVmmetI> PlanC ° f eXtemaJ ta " k grid W ' h a P 0 *"''* 1 in ^nt unci a flat axis at 

wh r e A >' " V ,+ l “ ***•' “ the user - s P e cified arc length for marching 
and is the surface area element. The radial point distribution 

no^Tt th , £ T 1Cngth betWCen P° ints in the Erection 

normal to the body surface. In the most general case, s is a Junction of L h 
and / since each point on the body surface may be stretched to a different 
outer boundary location. (For example, in the case of a bodv at a positive 
angle of attack m hypersonic flow, one may wish to have the outer boundarv 
of the gnd farther away on the top surface than on the lower surface). 
Typically, the points are stretched awav from the bodv exponentially. If crrid 
point spacing control is required at the outer boundary as well as at the bodv 
surface, or if more uniform grid point spacing is required awav from the inner 
body surface, hyperbolic tangent stretching can be used. The ability to 
control the gnd spacing at each end of the domain is useful when multiple 
zones consisting of stretched points or uniformly-spaced points are desired 
Also, the presence of more points in the far field of a component <md helps 
to improve gnd-overlap regions for the chimera overset-grid scheme [9] ‘ 

In this land of volume control specification, if an initial distribution of 
points is highly clustered in f or p. then these points tend to remain highly 
clustered even far away from the body. In order to obtain a more uniform 
far- field distnourion, the volumes specified from (5.1) are averaged in € and 
P with each step taken in £. For example, the averaged volum^ AV , , can 
be computed as 

*•) + j 

X ( A W/ + + &Vj. k *. U + AV ,.*-u). (5.2) 
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where this is applied one or more times with each step in £ . A typical value of 
v Q that has been employed is 0.16. 

An alternative method for specifying mesh cell volumes is to form a grid 
about a ‘similar* but simple reference body for which the grid can be 
generated analytically, and to use appropriately scaled cell volumes from this 
reference grid for the more complex problem. This method is described in 
more detail in [13]. 


6. SMOOTHING 

Since central differencing is used in the numerical scheme, artificial 
dissipation terms are added to control oscillations arising from the odd-even 
uncoupling of grid points. An equally important function of the added 
dissipation is to control the smoothness of the resulting grid. The form and 
magnitude of the added dissipation are very important in shaping the grid 
quality. 

The parameters 6 ^ and 6 V on the left-hand side of the grid generation 
Equation (3.4) can be thought of as a type of smoothing in the f marching 
direction. In the f and 77 directions, it is adequate for the purpose of 
controlling the smoothness of the grid to use just second-order smoothing. 
The form of this smoothing is described below. 

For simple surface topologies, such as an ellipsoid that has a convex profile 
in all directions, a constant dissipation coefficient is sufficient to provide good 
grid quality. However, body surfaces encountered in aerodynamic applica- 
tions are frequently much more complex with combinations of sharp convex 
and concave comers. Hence, the dissipation, which tends to reduce grid 
curvature and orthogonality, must be applied more selectively. The reduction 
of grid curvature can be advantageous in concave regions and detrimental in 
convex regions. In order to prevent grid lines from converging and crossing, 
the dissipation has to be relatively high in concave regions (e.g., in the region 
between the Orbiter fuselage and the wing root in Figure 3). Conversely, the 
dissipation must be kept small near the body surface and sharp convex 
comers (such as the convex comers of the IEA box in Figure 4), otherwise 
the resulting grid spacing in the £ direction would become too reduced or 
even negative. A spatially uniform dissipation coefficient is unable to satisfy 
all of the above requirements. Spatially-varying dissipation coefficients have 
been used previously, but they have only accounted for variations of mesh 
size, specifically, C”*A and C~ l B variations. A spatially- varying form of the 
dissipation coefficient that has worked very well for many cases is described 
below. 

The explicit second-order dissipation D t added to the right hand side of 
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Fic. 3. Plane on the back section of the Orbiter. 


the equations is given bv 

D ' " ( 6 . 1 ) 

with 



Fic. 4. Plane through IEA box on the ring of the solid rocket booster. 
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where e c is a user-specified constant of 0(1), and are approximations 
to the matrix norms ||C" l All and ||C~ 1 BIL respectively, given by 


N f 


• , 2 , 2 

f j + ft + =f 

*! + y't + z ( 


N. 


ft* ~ z i 


*i + w + =; 


and Re , R are the dissipation coefficients given by 


R( - S, «, = ^l d j.k.l a j.k.l‘ 


(6.3) 


(6.4) 


The dissipation coefficients consist of three functions that provide different 
wavs to automatically adjust the local dissipation appropriately, depending on 
the local grid topology: 

(1) A scaling function S h which varies with normal distance from the body 
surface. 

(2) A grid point distribution sensor function, dj k , or dj k h depending on 
the direction, which senses mesh convergence cased on the distances be- 
tween neighboring grid points. 

(3) A grid angle function, af kl or aj kl , depending on the direction, 
which senses mesh convergence based on the angles between neighboring 
grid points. 

The form of the scaling function S, is given bv 



(/-l) 

(L,. - i) 

(^n, - 1) 

(L», - 1) 


2 * / < hnns • 

K™, + i*i <i, 


(6.5) 


where l^, is the number of points in the l direction and / trjm is restricted to 
the range (f, 1] X With df kl and d] kJ defined by (6.Sa. b). is set 
to l when one or both of the following is true: 


max 

j- k 


d f.k.t ~ d f.k.i-\ 

j. k 


max 

J k 


clj.k.t - m f d j.k.t-\ 

J • " 


< 0. 

< 0. 


(6.6a) 

(6.6b) 


Once Z lran , is located, the above tests are not performed for / > l. 
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The purpose of the scaling fraction S, is to guarantee small dissipation, 
and hence, grid orthogonality near the body surface. As one moves away from 
the body surface, dissipation is increased since grid lines may begin to 
converge in concave regions and some loss of grid curvature in convex regons 
is no longer a problem. Away from the body surface, the grid point distribu- 
tion sensor function and the grid angle function alone are sufficient to 
provide the appropriate amount of dissipation. Hence, the influence of S, is 
removed by setting it to a constant at some location l traM away from the body. 
It is found that a good criterion for locating is when the convergence of 
local grid lines is slowing down in some sense given by (6.6a, b). For all the 
cases encountered so far, it is sufficient to test for for l > For 

some cases, ?t™ ns may be less than I max , which reduces loss of orthogonality 
near the outer boundary; while for other cases, may have to be equal to 
^moi » w hen increasing values of dissipation are needed all the way to the outer 
boundary to prevent grid lines from convergng. 

The forms of the gnd point distribution sensor functions , , d v , , are 
gvenby J J ' 

jf.k.i = max [(^/*./)' /s \o.i , 

( 6 . 7 ) 

jj.k.t = max [(^7.t./) 2 /s '.o.i], 

where 


J£ _ I* j+ l.k.l— 1 r j.k.l-l 

d J kl ^ ~ 




' r j~l.kJ * r j.kj\ + 10-1.*-# ~ F j.k. 


(6.8a) 


Jv _ r j.k.l-l\ + I*;.*-!./-! r j.k.l - l 

a j. k . / ^ ^ — 

“ 7 j.k.i\ + \ F j.k-u - r J kJ \ 


(6.8b) 


The distribution of grid points in the £ and 77 directions are monitored by the 
functions df k t and dj k h respectively. The quantity df k , is the ratio of the 
distances between a gnd point and its neighbors in the £ direction at level 
(/ 1 ) to that at level /. This ratio is high in concave regons and hence, 

more dissipation is provided here. It is of order one or smaller in flat or 
convex regons, where less dissipation is needed. Similarly, dj k , represents 
the corresponding quantity in the 77 direction. The quantities df k ; and df k , 
are constructed from df k , and df kh respectively, which are*' raised to the 
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power 2/S; in order to counteract, the small value of S ( close to the body 
surface. Also, the values of df k 7 and dj ki are limited from becoming too 
low by a limiter of 0.1. Note that a grid point distribution sensor function 
based on cell area ratios is not as effective since the grid lines could be 
converging in one direction but not the other and the cell areas would not 
decrease very much. 

The grid angle functions af kh aj k{ are more conveniently defined in 
terms of the following unit vectors. Let the vectors pointing in the plus and 
minus £ directions at grid point (j y kj) be represented by ry and r*~ , 
respectively, where 


0 


r j* 


l.k.t 


r j,k.b 


r j-\.kJ r j,k.b 


(6.9) 


and let rj and rj be the respective unit vectors for r* and r ~ . Similar 
expressions are defined for r k * , r k “ and r k , r k for vectors and unit vectors 
pointing in the plus and minus 77 directions, respectively at grid point 
( j % Jk, l). The local unit normal hj k l based on the cross product of the above 
unit vectors is given by 




(£ - £) x (K - £) 
K? - K) * - f *')l ' 


(6.10) 


The cosine of the local half angle oc j k j in the £ direction is then given by 

«> s a j.k.i = "j.k.i ■ 0* = n j.k.i ■ ■ ( 6n ) 

The grid angle function aj k ( is then defined as 

l/(l - cos 2 ot Jkl ) if 0 , 

(612) 

1 — < a j.k.i < "■ 

A similar expression is defined for the grid angle function oJ kl in terms of 
the cosine of the local half angle k j in the rj direction. 

The smoothing provided by the grid point distribution functions in (6.7) 
has to be modified locally at grid points located at very sharp concave 
comers. These comers are detected by computing the half angles a and (3 
subtended by neighboring grid points in the £ and 17 directions respectively 
(see (6.11)). At a severe concave comer point, extra dissipation is required to 
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-I 

prevent crossing of grid lines. The functions of u and aj k , serve to provide 
the appropriate local modifications to the dissipation required at these comer 
points. For a typi^l mesh, the values of af kJ and af kJ are close to unity at 
most gnd pornts. Figure 5 shows the grid marching out of a concave an<de of 
20 aegrees. Note that orthogonality is still maintained for the first point off 
the surface. Grids for concave angles down to 5 degrees have been obtained 
with the above algorithm. 

A minor local refinement at / = 2 of the above dissipation scheme for 
sharp comers is given in the Appendix. 


7. METRIC CORRECTION 


Often a body surface will have a comer or edge region that is a discontinu- 
ity to the hvperboUc grid equations. Nevertheless, the spatiallv-varving dissi- 
pation described m Section 6 will typically generate a sufficiently smooth <nid 
at either convex or concave comers if the grid spacing to each side of°the 
comer is approximately equal. If the surface grid possesses comers with 
uneven grid spacings, dissipation alone is usually not sufficient to provide a 
smooth volume grid and additional remedies are needed. The metric correc- 
tion procedure described below generally provides a satisfactory treatment of 
comer discontinuities in all but extremely convex cases. 

With ^-derivatives as defined in (3.5).' the direction in which the <md will 
emanate from the comer is such that it is perpendicular to the line* joining 
the two neighbor points of the comer. For a comer with unequally spaced 
points this is clearly not a desirable feature (see Figure 6 a). In order to <mide 
the gnd out in a direction that bisects the angles at a point subtended bv its 
neighbors in both the £ and 77 directions, the derivatives x, u. - at the 
point are modified to rj, y' (t z \ as follows: ‘ 


i/f 


AV’ 

Det(C') 



<Jv z ( 


X l9(J 


(7.1) 



Fic. 5. A concave comer at 20 decrees. 
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F 1C. 6. Comparison of treatment of concave comer with unequal grid spacing*, (a) Without 
metric correction, (b) with metric correction. 


with 


o«(c-) - *k>‘ * (*; % - «;*>*• 


where 




'o 


\ ,r ; 


(7.2a) 


|r > 


and 


( x v =i) r = ^ (i F n + !?*’!) 


1 r; 


\n 


7rn - 77^ • ("-2b) 


While the ^-derivatives should be computed bv (7.1) near the bodv 
surface, the original method of computing these quantities should be restored 
away from the body surface. This can be achieved smoothlv bv 

(*<•%. =<) r “(l - *,)(*?•! /?•:?/ +»,(*}. tf'l) 7 . (7.3) 

where = 2" 1 and x'j. j/^\ z'J are computed by (3.5). The form of (7.1) is 
the same as (3.5) except that the and 77 -derivatives are replaced bv the 
corresponding primed quantities. These modified and 77 -derivatives are 
constructed in such a wav that the neighboring points of the comer appear to 


198 


WILLIAM M. CHAN AND JOSEPH L. STECER 

be of equal distance from the comer. The result of applying this procedure to 
a concave comer is shown in Figure 6b. 


8. TREATMENT OF CONVEX CORNERS 


Two methods are presented here that provide extra robustness at convex 
comers. They both involve switching from solving the hyperbolic grid genera- 
tion Equation (3.4) to some other equations at the convex comer ,point. The 
first method is an implicit averaging scheme. Instead of solving the hyper- 
bolic equations at a convex comer, the following average equation is solved. 


where 


A 0.» = i( + as) a o.*- (8.1) 


H Ar /.* ~ n Ar r^- 

+ A o-,4 

(8-2a) 

M„ A o.i = i( A o.i-i 

+ A 0.*-i)- 

(8~b) 


In other words, the marching increment at the comer is the average of the 
marching increment of its four neighboring points. The form of the above 
scheme can be made compatible with the hyperbolic grid generation scheme 
given by (3.4) by approximate factorization. The equation to be solved at 
convex comer points is then 

(l-iM f )(l-j^)AF=0. (8.3) 

At l = 2, the switch to solve (8.3) is performed if the external angle of a 
convex comer in either the £ or t j direction exceeds 240 degrees. As the grid 
marches out in Z, the switch to return to the normal hvperbolic scheme is 
performed when the minimum of cos a i k and cos )3. k becomes larger than 
about —0.2. Although some factorization error is present, this implicit 
averaging scheme has worked well for a variety of convex comers, including 
some cases where the normal scheme has failed. One such case is illustrated 
by the NACA00I2 airfoil shown in Figures 7a and 7b. The point of this 
example is not to show that an O-grid should be used for the airfoil but 
simply to demonstrate the ability of the scheme to produce high quality grids 
around sharp convex comers. 

The implicit averaging scheme described above possesses two desirable 
properties. It can be shown geometrically that if the grid marches out' 
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(a) (b) 

Fit;. 7. Airfoil with sharp trailing edge and O-grid topology, (u) Far view, (b) close up view 
of sharp trailing edge region. 

orthogonally at the neighbors of the convex comer, then the distance marched 
out at the comer point is always smaller than the distance marched out at the 
neighboring points. Moreover, the distance marched out at the comer point 
becomes smaller as the convex comer becomes sharper. This is very useful in 
helping to bend the neighboring grid lines towards the sharp comer as the 
grid is marched out awav from the body surface (see Figure 7b). The second 
property is that the grid marching out from the comer point will bisect the 
angle at the comer provided the neighboring points march out the same 
distance and in symmetrical directions to each other. 

An alternative but potentially more robust method than the implicit 
averaging scheme above is described below. The exact location of the grid 
point in the next marching step out from a convex comer is predicted in 
advance. The predicted point is located by marching the grid a distance of 
&Sj k f in the direction given by the angle-bisecting unit normal defined bv 
(i^ , y' , z\) T in (7.1). The marching distance is scaled as 

A $ j k j = A Sj k j X min(sin a } k ,,sin (8.4) 

where A Sj k , is the user-specified arc length in the normal direction (see 
Section 5). The scaling with the sine of the half angle causes the comer point 
to march out a smaller distance than its neighbors, thus helping to bend the 
neighboring grid lines towards the comer. Thus, AF at the convex comer 
point can be computed in advance and combined with the grid generation 
Equation (3.4). From numerical trials with different geometries, it was found 
that this procedure also works well for many types of convex comers. 
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9. GRID QUALITY CHECKS 

Two grid quality checks are performed in the hyperbolic grid generation 
code. The first check is a cell volume check by decomposing each cell into six 
tetrahedra [7], This test is only passed if the volume of each tetrahedron in 
the cell is positive. The second check involves computing the Jacobian of each 
cell using a finite volume scheme employed by typical flow solvers. Negative 
volumes or Jacobians are reported to the user and smoothing parameters can 
then be adjusted appropriately to remove the bad cells. 

10. RESULTS AND APPLICATIONS 

The vectorized version of the current hyperbolic grid generator runs at 
142 megaflops on the CRAY-YMP and requires about 9.7 microseconds of 
CPU time per grid point (i.e., about a hundred thousand points per CPU 
second). As an example, generation of the largest grid in the Space Shuttle 
launch^vehicle grid system, the Orbiter grid for flight Revnolds number 
(98 X , , x 57 - 430122 points), takes 4.17 seconds of CPU time. This is 
about one to two orders of magnitude faster than typical elliptic <nid 
generators. " ' ° 

The various external components of the integrated Space Shuttle vehicle 
contain a wide variety of different geometric features that are found in many 
other applications. Hence, these geometric components provide good tests of 
robustness for the hyperbolic grid generation scheme described above. Some 
examples are given below. 

The following two examples show how the spatiallv-variable dissipation 
coefficient works at different types of comers. The first example shows the 
grid around a sequence of sharp convex and concave comers that appears at 
the IEA box on the attach ring of the solid rocket booster as in Figure 4. Low 
dissipation values are needed at the body surface and above °the convex 
comers to maintain orthogonality while high dissipation values are needed in 
the concave regions to provide grid smoothness. Figure 3 shows the smooth- 
ness of the grid in a large concave region, the wing root region at the back of 
the Orbiter. 

The robustness of the axis logic is demonstrated in the next example where 
the geometry is noncircular in the circumferential direction around the axis 
and that the distribution of grid points in this direction is nonuniform. This 
type of axis condition is present at the tip of the vertical tail section of the 
Orbiter as in Figure 8. The mixed extrapolation scheme with volume scaling 
is able to produce a smooth grid for this case. 

When the chimera overset-grid method is used on two grids whose body 
surfaces intersect each other, the grid points in the region around the 
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FlC. 8. Vertical tail grid showing parr of the surface grid near the tip. the 3D grid on the 
symmetry plane and a section normal to the symmetry plane. The axis is marked bv the thick 
dark line. 


intersection line which are common to both grids are left with no interpola- 
tion stencils. In order to remedy this problem, a collar grid can be introduced 
which covers the region around the intersection line [11], The most challeng- 
ing example tested by the grid generator so far is the collar grid, which covers 
the intersection region between the vertical tail and the Orbiter as in Figures 
9a, b. The methods used to generate the collar grid surface are explained in 
[9]. The surface of the collar grid is made up of two parts. The top part lies on 


(a) 


(b) 



Fic. 9. Views of the collar grid surface, (a) Position relative to vertical tail and Orbiter. (b) 
close up view relative to vertical tail and Orbiter. 
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the surface of the tail down td'the intersection line with the Orbiter. The 
lower part has two sections. The first section starts at the intersection line 
with the tail and then follows the top surface of the Orbiter. The second 
section folds over the back of the Orbiter and follows the backward-facing 
aft-bulkhead of the Orbiter. The difficult feature of this geometry is the 
presence of a region where grid lines are concave in one direction and convex 
in the other. Slices of the 3D grid viewed from the front and back ends of the 
collar grid are shown in Figures 10a, b. 

The final example is taken from the telescope grid for the SOFIA vehicle 
[ll. The SOFIA is a modified Boeing 747 with a telescope mounted inside a 
cavity on the upper surface of the plane. The telescope is topologically similar 
to a hollow bowl with a truncated cylinder in the middle of the inside of the 
bowl. Figure 11a shows the surface geometry (shaded) for half the telescope 
and slices of the 3D grid The external surface of the telescope consists of 
both the outside and inside of the bowl together with the middle cylinder. 
The symmetry plane of the 3D grid is shown in Figure lib. These figures 
show that, although the surface grid may not possess the sufficient number of 
points to resolve the detailed flow structures, the grid generator is able to 
produce a smooth grid over the complex combination of concave and convex 
comers using the techniques described in Sections 6, 7 and 8. The outer 
boundary of the telescope grid need not be placed far away from the body 
surface since the entire telescope grid is surrounded by a larger cavity grid 
and communication between the two grids is achieved via the chimera 
overset-grid scheme. 
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FiC. 11;l Surface geometry and slices of the 3D grid for the SOFIA telescope. 



Fie. Ilb. Symmetry plane of the 3D grid for the SOFIA telescope. 


11. CONCLUSIONS 

A robust three-dimensional hyperbolic grid generation scheme that is able 
to produce high quality grids for a wide variety of geometries has been 
presented. The improved robustness and the speed advantage of the scheme 
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have made it extremely attractive for users of chimera flow-solvers. Since the 
scheme is fast, the user can readily adjust the input parameters to fine-tune 
the grid quality. 

The use of a spatially-varying dissipation coefficient based on distances and 
angles between neighboring grid points gives the grid generator the ability to 
cope with geometries that are more complex than before. The grid angle 
bisecting property is provided by the metric correction procedure near the 
body surface. Sharp convex comers are automatically detected and the grid 
generation equations are altered at these comers to further enhance smooth- 
ness and robustness. 


APPENDIX 

In order to guarantee orthogonality near the body surface, the dissipation 
coefficients R { and R v described in Section 6 are made to be zero every- 
where at / = 2 through the scaling function S,. However, near concave arid 
convex comers, some dissipation has to be restored to maintain smoothness. 
This can be accomplished by introducing the blanking function b- k , which 
multiplies and R v at / = 2. The blanking function is zero everywhere 
except near concave ( a ; k or Rj k < it/ 3) and convex ( a, k or R. k > 2-/3) 
comers. For a concave or convex comer in the ^-direction, we set 

b j.k = K . b JtUk = 0.5, bj. ik = 0.25, (A .1) 

while for a concave or convex comer in the rj-direction, we set 

b j.k = K, b Jkj:l = 0.5, b j.k± 2 = 0.25, (A. 2) 

where b e = 1 for a concave comer and b c — 0 for a convex comer. At convex 
comers, dissipation at the neighboring points to the comer is still needed, but 
the dissipation at the comer point itself should be set to zero in order to 
produce the desired effects for the special schemes described in Section 8. In 
constructing bj k , it is assumed that successive comers in a coordinate 
direction are separated by at least four points. This is a reasonable assump- 
tion if one wishes to provide sufficient resolution for the flow around such 
comers. 
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ABSTRACT 

The application of CFD techniques to large 
problems has dictated the need for large team efforts. 
This paper offers an opportunity to examine the 
motivations, goals, needs, problems, as well as the 
methods, tools, and constraints that defined NASA’s 
development of a 111 grid/16 million point grid 
system model for the Space Shuttle Launch Vehicle. 
The Chimera approach used for domain 
decomposition encouraged separation of the complex 
geometry into several major components each of which 
was modeled by an autonomous team. ICEM-CFD, a 
CAD based grid generation package, simplified the 
geometry and grid topology definition by providing 
mature CAD tools and patch independent meshing. 
The resulting grid system has, on average, a four inch 
resolution along the surface. 
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1. INTRODUCTION 

This paper emphasizes the culmination of many 
separate grid generation efforts that, together, 
provided a grid system for the Space Shuttle Launch 
Vehicle (SSLV) that has a very high degree of 
geometric fidelity. The distributed nature and large 
scale of this project provided special challenges that 
are likely to become increasingly typical of 
Computational Fluid Dynamics (CFD) analysis. 

The SSLV has an extremely complex geometry for 
an aerodynamic vehicle. The four main SSLV 
components, the Orbiter, the External Tank (ET), and 
two Solid Rocket Boosters (SRBs), are sketched in 
Figure 1. The proximity of these bodies makes 
accurate modeling of the flow field and aerodynamic 
loads extremely dependent on the accuracy of the 
geometric model. 

Initial efforts to model the flow field of the SSLV 
began with a coarse grid that severely simplified the 
geometry of the flight vehicle 1 ' 2 . The importance of 
accurate geometry representation became readily 
apparent as geometry modifications improved 
comparisons between simulation results and wind 
tunnel data 3 ' 5 . For example, one relatively small 
modification that added the integrated electronics 
assembly and attach ring to the SRB caused the 
simulated flow under the Orbiter wing to change 
enough to increase the wing-load by approximately 
twenty-seven percent of the maximum allowable 
load 3 . Several other components of the flight vehicle 
geometry, such as the ET/SRB and ET/Orbiter attach 
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hardware, also block the flow between the Orbiter, 
ET, and SRBs. Unfortunately, the crudeness of 
previous grid generation tools and lack of concise and 
accurate geometric definitions hampered the 
refinement of the "preliminary" grid system. Much of 
the preliminary geometric model had been defined 
with available drawings and photographs that were 
often lacking in detail. These approximate surface 
definitions had been encoded and discretized within 
grid-specific FORTRAN programs. Typically, the 
resulting grids were then exposed to a tedium of small 
perturbations and "smoothing" operations as a final 
grid was evolved. This "coding and evolution" process 
was too time consuming, cumbersome and laborious to 
be considered for major refinement of the preliminary 
model. 




Figure 1. Front, Rear, Top, and Side view drawings 
of the SSLV configuration. 

The development of the present enhanced- 
resolution SSLV grid system began with a study of 
geometric modeling and grid generation requirements. 
Evaluation of available tools and data showed the 
desirability of a grid generation tool that could 
utilize the mature capability of Computer Aided 
Design (CAD) techniques and geometric databases. 
Several subcomponents of the launch vehicle had 
already been modeled, with many man-years of 
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effort, by other groups and were available. 
Furthermore, these CAD models represent the most 
accurate and complete geometric descriptions 
available for the launch vehicle. Several different 
CAD-to-grid generation tools were evaluated but most 
had only limited abilities for generating ordered 
meshes. Control Data Systems' software package, 
ICEM-CFD 6 ' 7 , provided the only CAD interface for 
grid generation that allowed an ordered mesh to be 
generated independent of the orientation of the 
underlying geometric patches. 

The grid system was generated for use with the 
flow solver OVERFLOW® that utilizes the Chimera 
domain decomposition approach 9-12 . The grids for 
each geometric component of the multi-body geometry 
were generated independently and overlapped to 
produce a complete grid system. During flow 
simulation the boundary values of each overlapping 
grid are updated through interpolation. Only the 
local resolution and amount of overlap between grids 
need be considered when interfacing two adjacent 
grids. Grid generation is not limited by the point-for- 
point or cell-face matching associated with "Block 
Zone" methods. This allowed a high degree of 
flexibility in determining the resolution and topology 
appropriate for modeling each component of the 
geometry. 

As with many other large-scale CFD simulation 
efforts, the Space Shuttle Launch Vehicle ascent 
simulation was accomplished in a team effort that 
was spread over several locations and organizations. 
Although full group discussions were conducted 
periodically, it was virtually impossible for each 
member to remain cognizant of all decisions. It proved 
vital to have a method of breaking the problem into 
separate tasks with limited interfacing requirements 
so that each group was encouraged to progress 
autonomously whenever possible. 

The approach, criteria, tools, and limitations 
associated with the enhanced-resolution SSLV grid 
system are discussed below in the geometric definition 
and grid generation sections. 

2. GEOMETRY DEFINITION 

The search for data that defined the outer mold 
line surface geometry of the Space Shuttle Launch 
Vehicle revealed several sources including: Space 
Shuttle Master Dimension Specification 13 , Interface 
Control Document # ICD-2-00001 14 , numerous 
photographs, literally thousands of drawings and 
several CAD and numerical descriptions. 
Responsibility for subcomponents and related data 
was scattered over several NASA centers and 
contractors. Although this information was 
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available for use by NASA, creating a concise and 
appropriately accurate model posed to be a major 
undertaking. The level of detail varied from source to 
source. Manufacturer's production drawings 
represented the most accurate definition of the actual 
vehicle components but were far too numerous to 
compile for an entire vehicle in the time allotted. 
The Master Dimension Specification contains a 
detailed description of the Orbiter but left the ET, 
SRB, and protuberance definitions to be described by 
ICD-2-00001. Both documents were encumbered with 
a complicated revision process where a current surface 
geometry could only be determined by superimposing 
all previous revisions. However, these documents did 
provide accuracy appropriate for spot checking and 
verifying other models. The decision to use a CAD 
based grid generation package provided a medium for 
compiling definitions and a means of utilizing existing 
databases that were previously inaccessible. CAD 
information was collected from several institutions 
and compiled into major component CAD models 
representing the Orbiter, External Tank, and Solid 
Rocket Boosters. These three major component models 
were verified, and supplemented where necessary 
with information from drawings, photographs, and 
flight vehicle inspections. 



Figure 2. Front view solid shaded image of the 
enhanced resolution SSLV surface grid. Note the 
resolution of the cable trays and feed line on the 
ET, the circumferential plant and field joints and 
cable tray on the SRB, and the Orbiter canopy and 
forward bipod. 

Each of these major components was developed 
independently with only occasional updates between 
groups. Geometry definitions of the separate 
protuberances and components of the ET were created 
by Martin Marietta at Michoud. The highly detailed 
CAD descriptions of the ET were modified to remove 
internal features, converted to surface models, and 
composed into an outer moldline surface model by 
Lockheed Engineering and Sciences Company at 
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Houston, Texas. CAD definitions of portions of the 
Space Shuttle Orbiter geometry were created by 
Rockwell International Space Systems Division at 
Downey, California and at Palmdale California. 
Composition of the Orbiter models into a CAD model 
suitable for CFD was performed by Rockwell-Downey 
and NASA Ames Research Center. The SRB geometry 
definition was created and composed with 
cooperation amongst Rockwell-Downey, NASA Ames, 
and NASA Johnson Spacecraft Center. The following 
is a detailed discussion of the efforts associated with 
each major component. 

2.1 External Tank 

The true surface geometry of the ET includes many 
protuberances such as feedlines, cable trays, fairings, 
and large structural members (for attaching the ET to 
the Orbiter and SRBs). Each of these protuberances is 
exposed to the exterior flow field during launch, and 
many present large cross sections normal to the flow 
field. The enhanced-resolution flow field simulation 
needed to be accurate enough to compare well with 
flight load and flow field measurements. A target 
surface resolution of four inches for the overall launch 
vehicle model was determined largely by the 
resolution required to capture the effects of these ET 
protuberances. 

Defining the geometry of all tank protuberances 
greater than four inches using the tools and data 
available for the preliminary grid system would 
have been an insurmountable task. The complexity of 
the geometry, the continuous revision process, and the 
accuracy required made it difficult to identify and 
locate accurate and current information. This was 
exacerbated by the difficulty in defining the 
nonuniform layer of protective thermal foam over 
each component of the ET. Some outer moldline 
definitions were unavailable due to the operational 
approach of defining the thickness of the foam with 
large tolerances and allowing the line worker to 
define the final shape. 

Fortunately, Martin Marietta at Michoud had 
already undertaken an effort to update the definition 
of each protuberance of the ET. The resulting models 
were extremely detailed three-dimensional 
wireframe CAD definitions. Each outer moldline 
definition included descriptions of the underlying 
metal structure and each layer of thermal protection 
foam that is applied during assembly of the ET. 

A formal request was made for the electronic 
transfer of these protuberance models and a 
description of the tank surface itself from Martin 
Marietta to NASA JSC. The transfer required a 
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conversion from the Computer Vision CAD system, at 
Martin Marietta, to the Control Data System ICEM- 
CFD CAD system, at NASA-JSC. This was 
accomplished using IGES format translators within 
each code (see section 2.4). To expedite the transfer, 
"clean up" of the models (to remove internal details 
and small elements such as bolts and lockwire that 
were obviously below the anticipated resolution of 
the computational grid) was limited to those features 
that were easily removed prior to transfer. 






Figure 3. A rear view comparison of three models 
and the SSLV ET/Orbiter aft attach structure. The 
subfigures are: a) a solid shaded image of the 
preliminary grid system, b) a photograph of the 
IA-613 wind tunnel model, c) a solid shaded image 
of the enhanced resolution grid system, and d) a 
photograph of the actual flight hardware. 
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As the tank subcomponent models were received, 
the wireframe definitions were used to create an outer 
moldline surface model of the ET and the ET/Orbiter 
attach hardware. In the early stage of the 
redefinition, the surface geometry was not modified 
beyond the level convenient for creating the 
component surfaces expediently. The outer mold line 
surface geometry was generally left as it was 
originally defined in the Martin Marietta CAD 
definition. Considerable effort was required to get all 
components into a common coordinate system and to 
cycle through the inevitable inconsistencies 
encountered when assembling any large collection of 
parts. 

The intermediate result of collecting these models 
was a very detailed CAD surface model. A 
considerable amount of surface approximation was 
required to remove fine surface details and simplify 
the topology to allow discretization with a 
reasonable number of points while maintaining the 
targeted four inch surface resolution. However, the 
process of eliminating unwanted detail resulted in a 
different approximation than would have been 
achieved if an approximating surface had been 
created from scratch. It was easier to retain detail 
and the accurate surfaces than to create a 
geometrically simpler approximating system. Thus 
the resulting grid surface is highly faithful to the 
true geometry except where the topology was severe 
enough to warrant the effort required to change it. 

The cable trays on the ET surface are modeled as a 
single constant-height ridge in the tank surface that 
ignores the many separate cable tray and repress line 
tie down brackets. The L02 cable tray (that runs from 
the intertank region, below the Orbiter nose, to a 
fairing at the forward tip of the ET) can be seen in 
Figure 2. The multi-faceted intertank region of the ET 
is modeled as a smooth cylindrical surface with 
circumferential ramps modeling the transition from 
the intertank region to the rest of the tank surface. 
The L02 feed line, tie down brackets, and forward 
fairing is modeled as a bump in the ET surface that 
can be seen in Figure 2 beginning below the Orbiter 
nose and running back to the aft attach hardware. 
Figure 3 shows a progression of fluid dynamic models 
of the ET/Orbiter aft attach hardware model in 
comparison with the flight hardware. The 
preliminary (symmetric) grid system, shown in Figure 
3a, captures most of the large elements that block the 
flow but does not include the intersections between 
bodies. The wind tunnel model. Figure 3b, includes the 
diagonal strut and umbilical cable trays but 
underestimates the size of each component by the 
thickness of the thermal protection foam. The 
enhanced resolution grid system. Figure 3c, can be seen 
to be the most accurate model of the actual flight aft 
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attach hardware. Figure 3d. Comparing the LH2 feed 
line region model on the left side of Figure 3c to the 
flight hardware in Figure 3d shows that virtually all 
the features are represented. However, the L02 
region, on the right side of Figure 3c, is missing some 
detail. The difference between these two models 
reflects that on the LH2 side each component was 
relatively exposed while many of the features on the 
L02 side were so tightly clustered that many were 
sheltered from the oncoming flow. As a result, some 
components on the L02 side were deleted from the 
model and several were combined to form one L02 
"blob" by merely filling in the small remaining gaps. 
In effect, the approximation fits the actual surface 
like shrink-wrap. 

2.2 Obiter 

The Space Shuttle was designed in the early 
seventies using parametric equations to define the 
cross sections of the Orbiter outer mold line at discrete 
locations along the fuselage. The Numerical Design 
group at Rockwell-Downey maintains this geometric 
database via the Orbiter Master Dimension 
Specification Program. The output of this program 
yields point data that models a smooth surface 
description of the Orbiter that is devoid of die gaps 
and depressed regions where riles have been replaced 
by thinner thermal blankets. Also missing are the 
variations found in each of the Orbiter vehicles due 
to manufacturing tolerances. Both the preliminary 
and enhanced grid surface models of the Orbiter were 
initially based on surface points supplied by 
Rockwell's Numerical Design Group. 

The preliminary grid system was developed from 
surface points that defined the Orbiter cross section at 
107 specific stations, or approximately every ten 
inches, in the streamwise direction. A FORTRAN 
program was written to create spline curves at each 
streamwise station and to interpolate between the 
splines to create surface grid points in a structured 
mesh. This model was accurate but lacked detail 
towards the rear of the Orbiter and the surface grid 
topology was difficult to manipulate. 

The initial enhanced-resolution CAD geometry 
was developed at Rockwell-Downey, using ICEM- 
CFD, from a much larger set of surface points. Point 
data was obtained at each half inch interval along 
most of the fuselage and at each one hundredth inch 
interval near regions where the geometry changed 
rapidly such as the nose, canopy, Orbiter 
Maneuvering System (OMS) pods, and elevon regions. 
B-splines were created from select point data at each 
cross section and B-spline surfaces were created from 
select curves. Sixteenth order surfaces and tight 
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tolerances were used to insure accuracy. Sharp comers 
were modeled with multiple segment B-splines. This 
tedious interactive process has now been automated 
within ICEM-CFD. The resulting model was smooth 
and fairly free of surface undulations. 

However, the need for accurate definitions of the 
elevon gaps and wing tips led to a search for even 
more accurate models of these regions. An effort was 
underway to create a CAD description of portions of 
the Orbiter surface at NASA Kennedy Spacecraft 
Center and Rockwell-Palmdale to support tile 
fabrication activities. A formal request was made to 
obtain the CATIA CAD system models in IGES 
format. These detailed models were used to improve 
the CAD definition of the Orbiter lower wing 
surfaces, elevons, elevon flipper doors, elevon gaps, 
and wing tips. 



Figure 4. Rear view solid shaded image of the 
enhanced resolution SSLV surface grid. Note the 
resolution of the elevons and flipper doors on the 
Orbiter wings and the many concentric plant and 
field joints near the base of the SRBs. 

Orbiter base components such as the Space Shuttle 
Main Engines (SSMEs) and the body flap were created 
by various members of the team and combined into a 
consistent model at NASA Ames and Rockwell- 
Downey. Cooperation between these sites was greatly 
enhanced by the use of a common CAD package. The 
resulting enhanced grid system surface for the Orbiter 
is presented in Figures 2,4, and 5. The smooth forward 
portion of the Orbiter and canopy can be seen in Figure 
2. Figure 4 shows the complexity of the rear of the 
Orbiter wings and base. The wing models include 
inboard and outboard elevons, gaps, and upper surface 
flipper doors. Figure 5 provides a close up view of the 
Orbiter base. Note that the SSMEs, OMS pods, and 
Reaction Control System (RCS) extensions are all 
included. The OMS engine bells were excluded since 
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they do not fire during ascent and are located in the 
separated region caused by the OMS pods. 



Figure 5. Close-in rear view solid shaded image of 
the enhanced resolution SSLV surface grid. Note 
the detailed model of the Orbiter base region and 
aft ET/Orbiter attach region. 

2.3 Solid Rocket Booster 

No detailed CAD descriptions of the Solid Rocket 
Boosters and ET/SRB attach hardware similar to 
those available for the External Tank and Orbiter 
were found. Instead, the surface models were created 
from analytic definitions derived from drawings 
published in ICD-2-00001. However, most large 
protuberances on the SRB surface are included: The 
segmentation rings (plant and field joints) can be seen 
in each of Figures 2,4 and 5; The external cable tray 
can be seen in Figure 2; while the integrated 
electronics assembly box can be seen best in Figure 5. 
The current grid system models the Redesigned Solid 
Rocket Booster (RSRB) and work is in progress to 
model the Advanced Solid Rocket Booster (ASRB). 


2.4 IGES Interface 

The Initial Graphics Exchange Specification 
(IGES) provided a standard interface between CAD 
systems, that was reliable most of the time, from 
several different manufacturers. The transfers of wire 
frame CAD models encountered very few problems as 
points, lines, arcs, curves, etc., were translated to and 
from IGES format. However, increasing numbers of 
surfaces, which have not been established within 
IGES, were included in later models. Trimmed 
surfaces and surfaces of revolution were particularly 
difficult to transfer consistently. Extremely high 
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order surfaces also caused problems when the 
capability of the receiving system was exceeded. In 
most cases, the offending surfaces were merely deleted 
by the ICEM-CFD translator and required recreation. 
Close contact with the vendor and Control Data 
Systems personnel provided work-arounds and custom 
fixes for specific problems. 

3. GRID GENERATION 

The capabilities of ICEM-CFD and Chimera 
domain decomposition combined to greatly improve 
the standard approach to grid generation. In this 
project, ICEM-CFD was used primarily to develop 
surface geometry and surface grids. User interaction 
with the CAD interface of ICEM-CFD greatly 
reduced the cycle time in the iterative process of 
defining and refining entities that control the 
geometry and topology of the grid. Volume grids were 
generated with the hyperbolic grid generation 
program HYPGEN. The three-dimensional flow field 
domain was decomposed into overlapping regions in a 
Chimera approach that allowed each geometric 
component of the launch vehicle to be modeled in a 
separate effort. The program PEGSUS was used to 
interface the resulting overlapped volume grids. 
HYPGEN and PEGSUS are described in Sections 3.2 
and 3.3, respectively. 

g.l Surface Grid Dev e lo pm ent 

One key advantage of ICEM-CFD is its ability to 
develop a grid without being constrained by the 
orientation or extent of the underlying network of 
geometric patches. This "patch independent 
meshing" decouples the process of defining the 
geometry from the process of defining the orientation 
and distribution of the surface mesh. The creation of 
the CAD definition can proceed in a manner 
convenient for modeling geometric features while the 
topology and density of the mesh can be overlaid 
according to the needs of flow field features. 

Once the CAD surface geometry is complete a 
computational mesh can be defined. The interactive 
grid preprocessor module, MULCAD 7 , within ICEM- 
CFD uses a combination of base level CAD tools and 
menu prompts within a subprocess of the CAD 
environment to define the extent, topology, and 
clustering of the grid. The user has the convenience of 
alternating between geometric manipulation and the 
set up of grid control elements. This can be 
particularly useful in simplifying detailed geometric 
models to a level appropriate for the given discrete 
resolution and budget of grid points. 
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Figure 6. Qose up view of the surface grid of the SSLV showing the transition from an H-type to an O-type 
grid topology near the Orbiter nose. Also note the details captured within the ET surface grid. 


In this phase of generating an ordered mesh with 
MULCAD, a second network of geometric elements is 
created over the original geometry. The layout of 
this second network can be defined with the same 
CAD tools that were available for the geometry but 
oriented to explicitly control the grid topology. This 
approach complimented the capabilities of the 
Chimera domain decomposition scheme by allowing 
varied topologies to be overlaid on the same 
geometric model. Figure 6 shows an H-grid topology 
on the nose of the Orbiter smoothly transitioning to 
the O-type topology on the Orbiter body. 

The distribution of points along each control 
element is defined interactively and can be viewed for 
qualitative verification within MULCAD prior to 
computing the grid. The resulting hierarchy of grid 
control elements: grid control curves (called Edges) 
combine to define grid control surfaces (called Sub- 
faces) which combine to form volume or surface block 
definitions (called Domains). Typically, curves 
available from the preceding geometry definition 
phase are supplemented with curves created 
specifically for grid control then combined and given 


control attributes within MULCAD. The usual 
varieties of tools are available for distributing points 
along curves and surfaces, and within volumes. 

Experience with the preliminary grid system 
dictated which geometric components needed to be 
included in the enhanced-resolution grid system. The 
resolution required to include these components 
defined the four inch nominal grid resolution for the 
grid system. Of course, the grid was clustered tighter 
than nominal near regions with complex geometry 
and/or flow field features and relaxed near smooth 
geometry or low gradient regions. The final grid 
system has a surface resolution that typically varies 
from two to twelve inches. A minimum of three points 
were used to resolve geometric slope discontinuities, 
and five points were used where possible. Shocks are 
resolved by the flow solver with five grid points and 
central differencing with 2nd- and 4th-order 
smoothing. 

The large number of surface grid points required to 
capture the SSLV geometry is evident in the 
representation of the ET in Figures 6,7. In some cases, 
the number of points required to capture a geometric 
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Figure 7. SSLV surface grids in ET base region. Note how the domain decomposition simplifies the topology 
of each individual element and restricts the areas affected by local increases in resolution. 


element exceeded the number warranted for its 
inclusion. As a result, the geometry and/or topology 
was often simplified and the preceding processes were 
repeated. Components immediately downstream of 
larger components were assumed to be in separated 
flow regions that would cause them to have small 
effects on the flow field. Therefore the grid system 
most accurately depicts those that are directly 
exposed to oncoming flow. In general, any component 
smaller than a four inch diameter was not considered 
resolvable. 

Upon completion of the definition phase of 
surface mesh development, an input file is written 
from MULCAD. PADAMM 7 , the grid generation 
module for the ICEM-CFD package, is a batch 
program that uses the input parameters and grid 
control element geometry from MULCAD to generate 
an intermediate mesh and "project" it to the 
underlying CAD surfaces in a proprietary process. 
PADAMM provides both algebraic and elliptic tools 
for generation and smoothing of the projected surface 
grids. Experience showed that the resulting mesh 
faithfully followed the geometry of the CAD surface 


model but that special control was needed to maintain 
the intended distribution of points in regions of high 
curvature. When the control geometry oversimplified 
the underlying CAD surface model the resulting grid 
was not smoothly distributed. Additional control 
lines were added to simplify the projection by 
breaking the geometry into simpler pieces. However, 
neither the control line geometry nor the point 
distribution along the control line is adjusted by 
PADAMM during grid generation. Unfortunately, 
this effectively requires the user to guess the exact 
location of the controlling grid line within the mesh 
which can be a time consuming iterative process. 

Originally, PADAMM was executed as a batch 
process separate from the CAD environment. Towards 
the end of this project the ability was added to 
execute PADAMM from within ICEM-CFD. This 
reduced the cycle time for interacting with the 
geometry, topology, and distribution of the grid. The 
capability to generate volume grids within ICEM- 
CFD was not utilized in favor of the hyperbolic grid 
generator described in the section below. 
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Figure 8. Side view of the SSLV grid system with volume grids plotted on the pitch plane. Note that the 
surface oriented grids extend only a short distance before the medium resolution rectangular box grids form a 
transition to the coarse far field grid. The blank region behind the Orbiter is discretized by the SSME 
plume grids that do not lie in the pitch plane and were not plotted. 


3.2 Volume Grid Development 

The hyperbolic grid generator HYPGEN 1 s 
developed at NASA Ames Research Center was used 
to create the volume grids for several reasons. The 
hyperbolic marching scheme required one to two 
orders of magnitude less CPU time than typical 
elliptic methods. The user effort to set up the volume 
grid generation process was focussed on producing a 
high quality surface grid rather than being divided 
to define all six faces of the volume control block as 
with elliptic and algebraic schemes. The release of 
control over the extent and shape of the outer 
boundaries required by the hyperbolic method was 
appropriate for generating grids for the Chimera grid 
system. No constraints are imposed on the shape or 
location of grid outer boundaries beyond maintaining a 
one-cell minimum overlap between neighboring grids. 
On the other hand, tight control of clustering and 
orthogonality, which is provided by hyperbolic 
methods, was required for the grids to be used with 
the thin-layer Navier-Stokes approximation in the 
flow solver OVERFLOW. The enhanced robustness of 
HYPGEN allowed high quality grids to be produced 
with very few attempts for even highly irregular and 
nonconvex geometry. 


The user interface, UI 16 , provides a convenient 
display of HYPGEN control parameters on a Silicon 
Graphics workstation and speeds the development 
cycle time. Default values for the input parameters 
are automatically initialized and can be easily 
modified. UI can be run entirely locally on the 
workstation or the remote processing option can be 
invoked to allow the CPU intensive grid generation 
process to be run on a faster remote machine such as a 
CRAY. (This reduced the execution time for 
generating a volume grid from nearly an hour to a few 
seconds.) Viewing the resulting grid is provided by UI 
by creating a window running PLOT3D 1 ' and command 
files for setting the appropriate viewpoint and 
displaying default surfaces. 

Typical inputs for SSLV grids with no-slip 
surfaces were: a first cell normal spacing of .00016 
inches, an exponential distribution in the normal 
direction with a constant stretching rate in the range 
1.25 to 1.29, and a maximum normal extent of 20 to 60- 
inches. The very small normal spacing corresponds to 
y + =1.0 at ten percent aft on the launch vehicle with a 
flight Reynolds number of 2.7xl0 8 . A higher 
stretching rate and normal extent was used with very 
smooth geometry and somewhat lower values were 
required for stable operation of HYPGEN with the 
more irregular surface elements. 
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Figure 9. Front view of the SSLV grid system with 
volume grids plotted on a plane midway along the 
fuselage. Note again the short extent of the 
surface fitted grids. 

Figures 8 and 9 illustrate that the hyperbolically 
generated surface oriented grids can have topology 
and outer boundary shape independent from the 
surrounding domains. This made volume grid 
generation one of the quickest tasks. 

In some cases, a surface grid was incompatible 
with HYPGEN and a volume grid could not be 
generated properly. Typically, this could be easily 
resolved by simplifying the underlying surface or 
surface grid 

3.3 Domain Decomposition and Interfacing 

In order for a flow solver to update a solution on 
the interior of a grid each point on that grid's 
boundary must be provided with a boundary condition. 
In the case of those points that overlap other grids 
within the Chimera grid system interpolation is used 
to update the appropriate flow values. The grid 
interfacing program PEGSUS 18 , developed by 
CALSPAN at Arnold Engineering Development 
Center, preprocesses the grid system to determine the 
appropriate information to be provided to the flow 
solver for updating each grid boundary point of the 
system during simulation. An interpolation cell that 
surrounds the boundary point is located in a 
neighboring grid and bilinear interpolation constants 
are calculated based on the relative location of the 
boundary and interpolating cell points. The resulting 
"stencil" information is stored for use during flow 
simulation. PEGSUS was modified at NASA Ames to 
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provide for the demands of the viscous spacing and to 
accelerate processing time. 

The first step in dividing the flow field was to 
determine how to divide the complex surface 
geometry into topologically simple regions. Each 
region was made as large as possible to minimize the 
number of separate grids and interface overhead 
required to model the geometry. Larger grids also 
increase vector length and vector speed-up during flow 
solution. On the other hand, the 64 megaword in-core 
memory limit, and the 16mw optimum operational 
queue size of the Cray-Xmp computer at NASA JSC 
limited the number of points that should be included 
in each grid. 

Early estimates of the requirements of the flow 
solver OVERFLOW indicated that thirty-two words 
of memory were required per grid point. Thus the 
maximum volume grid size was limited to 500,000 
points. The tight initial cell spacing, cell stretching 
criteria, and normal extent target limited the 
maximum number of surface points that could be 
included in each grid to approximately 15,000. This 
number and the nominal four inch spacing criteria 
determined the surface area that could be modeled 
with one surface grid domain. 

The actual decomposition of a surface model into 
overlapping surface grids was further adjusted to 
reflect local detail and grid interfacing criteria. In 
the aft attach region grid sizes were limited to 
simplify the topology of each grid, see Figure 7. The 
struts, feedlines and cross members of the aft attach 
region were modeled as O-type grids that wrapped 
around the geometric element and flared out onto 
adjacent elements to provide the required overlap. 
This allowed grids to be produced with viscous 
clustering in only one index direction as is appropriate 
for the thin-layer Navier-Stokes flow solver. 

When two adjacent grids have similar resolution 
and topology where they overlap the Chimera 
scheme is merely a generalization of the Block Zone 
approach. This simple overlap can be easily 
constructed for simple geometry and provides a quick 
method of decomposing large smooth areas into 
appropriately sized grids with reliable interfaces. 
Higher local resolution can be utilized in one grid 
while maintaining similar topology at the interface. 
However, a change in resolution across interfaces is 
best limited to a factor of two especially near high 
gradient regions to avoid excessive interpolation 
errors in determining boundary values during the flow 
solution. 

It can be cumbersome to maintain simple overlaps 
between grids where the topology of the geometry is 
complicated. One key advantage of the Chimera grid 
scheme is that the topology of each grid is really 
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independent. Topologically simple grids can be used 
without regard to the surrounding system as long as 
the resolution at the interface between grids is 
comparable. This is exemplified by the rectangular 
box grids, that form the interface between the 
complicated surface oriented grids and the simple free 
field grid, see Figures 8-10. If grid points fall inside 
bodies or into regions where the resolution is not 
similar (and would not yield proper interpolation 
stencils) then the unwanted points are "blanked", 
flagged as not contributing, to the solution. The 
nearby points that are adjacent to the blanked region 
of the grid are treated as boundary points and 
updated through interpolation. The process of 
identifying points that should be blanked, "cutting 
holes", and the process of providing interpolation 
stencils can be accomplished using the program 
PEGSUS. 



Figure 10. Side view of the outer domain of the 
enhanced-resolution SSLV grid system. The blank 
region that extends to the exit of the far field grid 
is discretized by the SRB plume grids that do not 
lie in the pitch plane and were not plotted. 

The tools available within PEGSUS for cutting 
holes in a grid are based primarily on families of grid 
points from adjacent grids. Occasionally, "phantom" 
grids are created to cut holes that are not easily 
defined. These phantom grids are excluded from the 
resulting grid system since they do not contribute to 
the flow simulation. 

As a rule, each grid was generated such that its 
outer and hole boundary points overlap neighboring 
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grids by at least two grid cells. This insures that 
interpolation stencils would be based solely on points 
that were updated within the flow solver. Two 
overlapping stencils will continue to trade 
information that is independent of the surrounding 
flow simulation. A point that can be updated without 
relying on interpolated points is said to have a stencil 
with Quality, Q, equal to one. Quality values in the 
range 1>Q>0 are assigned based on the relative 
location of the point to the interpolated points within 
the stencil. The number of times that PEGSUS 
searches for stencils and the minimum quality that 
will be accepted during each search can be defined in 
the input. The 16 million point enhanced-resolution 
grid has 1.2 million points that require stencils. Most 
stencils are of quality equal to one while 2.2% and 
0.1% have qualities in the ranges 1>Q>0.6 and 
0.6>Q>0.2, respectively. Stencils with Quality less 
than .2 would exchange little or no new information at 
the overlapping boundary during iteration. Points 
that failed the minimum quality test are labeled as 
"orphans". The process of manipulating the PEGSUS 
inputs to satisfy recalcitrant orphans can be time 
consuming. 

In practice, it can be difficult to ensure good 
overlap with three-dimensional grids that have 
varying spacing and topology, especially near holes. 
If acceptable interpolation stencils could not be 
established, it was necessary to loop back to previous 
steps to modify the resolution, extent, and orientation 
of the grids. Test cases to gain experience in 
developing models that provided compatible 
interfaces helped to minimize the time spent 
reproducing the underlying grids. 

A serious impediment to the timely completion of 
this grid system was the lack of an interactive 
environment for the interfacing procedure. Roughly 
half of the grid system development time was 
consumed in this phase of the project. The relation 
between the three-dimensional hole boundary points 
and the surrounding grids was not always easily 
visualized. PLOT3D running on an IRIS graphics 
workstation was severely taxed even with its color 
display capabilities and fast local rotation. A three- 
dimensional viewing device would have been very 
helpful. 

CONCLUSIONS AND RECOMMENDATIONS 

An enhanced-resolution grid system, presented in 
top view in Figure 11, was developed to provide a 
highly accurate computational grid model of the 
Space Shuttle Launch Vehicle. It includes 111 volume 
grids in an overlapped grid system that contains 16 
million points. Each domain includes less than 
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500,000 points to allow Computational Fluid Dynamic 
simulations to be calculated on computers with 
moderate amounts of in-core memory. Most features 
larger than a four inch diameter have been included. 

The ability to utilize the CAD interface for 
creating and modifying the geometry and for access to 
existing geometric databases proved to be very useful. 
It would have been inconceivable to consider creating 
a grid system of this complexity using previously 
available tools. 

The Chimera grid interfacing scheme was useful 
in its flexibility both for decomposing the flow field 
geometry into separate domains and to divide the 
grid generation effort into separate autonomous tasks. 

Given the current trend towards more accurate and 
complete representations of complex flow fields, it is 
important to maintain and enhance the tools suitable 
for modeling and manipulating large systems of 
geometry, data, and grids. 
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Figure 11. Top view of the surface grid for the enhanced-resolution Space Shuttle Launch Vehicle, 
varied topology and density of the overlapped grids. 


Note the 
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